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Short range correlations are introduced using unitary correlation method in a relativistic ap- 
proach to the equation of state of the infinite nuclear matter in the framework of the Hartree-Fock 
approximation. The effect of the correlations in the ground state properties of the nuclear matter 
is studied. 



In this paper we introduce short range correlations in a relativistic approach to the description of nuclear matter 
for the first time. Although there are several procedures which may be used to introduce short range correlations into 
the model wave function, we work with the unitary operator method as proposed by Villars y). There are several 
advantages in using an unitary model operator. In particular, one automatically guarantees that the correlated state 
is normalized. The general idea of introducing short range correlations in systems with short range interactions exists 
for a long time @, H| but has not been pursued for the relativistic case. 

Nonrelativistic calculations based on realistic NN potentials predict equilibrium points that are not able to describe 
simultaneously the correct binding energy and saturation density; either the saturation density is correct but the 
binding energy is too small, or the correct binding energy is obtained at a too high density . In order to solve this 
problem a repulsive potential or density-dependent repulsive mechanism [l| is included. Due to Lorentz covariance 
and self-consistency, relativistic mean field theories |fj include automatically contributions which are equivalent to 
n-body repulsive potentials in non-relativistic approaches. 

In non-relativistic models the interaction arises from the interplay between a long range attraction and a very 
strong short range repulsion and it is indispensable to take short range correlations into account. In relativistic mean 
field models, the parameters are phenomenological, fitted to the saturation properties of nuclear matter. Short range 
correlation effects are may be included to some extent in the model parameters. However, we want to study the 
consequences of taking these effects into account explicitly • 

We start with the effective Hamiltonian as 

H = J ^ a (x)(-ia- V + /3M) Q/3 V/3(x) dx + ^ J vbl(x)^{y)V aPr 4\x - y\)yb s {y)^ p {x) dx dy (1) 

with 

V a/ 3 t7 s(r) = {(f) a p(P)i5V a {r) + {5 a p8 lS - a af3 ■ a lS ) V u (r) (2) 

where 

„2 p —m„r „2 p—m^r 

V„(r)=-f- , V u {r) = j L , (3) 
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and a are the Dirac-matrices. In the above, ip is the nucleon field interacting through the scalar and vector potentials. 
The equal time quantization condition for the nucleons is given by 

[tp a (x,t),ipp(^,t)^] + = 5 a pS(x-y), (4) 

where a and (3 refer to the spin indices. The field expansion for the nucleons if) at time t=0 given by [3] 

*W = ^ E [ U r$K,k + Vr{-k)c^_ J e**, (5) 

where U r and V r are 



For free spinor fields, we have cosx(/c) = M/e(k), sinx(fc) = |fc|/e(A;) with e(fc) = Vk 2 + M 2 . However, we will 
deal with interacting fields so that we take the ansatz cos\(k) = M*(k)/e*(k), sinx(fc) = \k*\/e*(k), with e*(fc) = 
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fc* + Af* 2 (fc), where fc* and M*(k) are the effective momentum and effective mass, respectively. The equal time 



anti-commutation conditions are 



C r,k> C lrJ+ = 5 rsh,k> = [ 5 r,fc> 5 !.p]+ • ( 7 ) 



The vacuum | 0) is defined through c r £ | 0) = c' -, | 0) = 0; one-particle states are written |fc, r) = - | 0); 
two-particle and three-particle uncorrelated states are written, respectively as \k,r;k\ r') = - c , | 0), and 
Ifc, r; fc', r'\ fc", r") = c 1 - - c , - I 0), and so on. 

1 r,k r' ,k' r" ,k" 

We now introduce the short range correlation through an unitary operator method. The correlated wave function 9] 
is | "J) = e tn |$) where |$) is a Slater determinant and is, in general, a n-body Hermitian operator, splitting into a 
2-body part, a 3-body part, etc.. The expectation value of H is 

= (gjgjg) = ($|e-'» H e in |<&) 

(*|*> ($|$> ' 1 J 

In the present calculation, we only take into account two-body correlations. Let us denote the two-body correlated 
wave function by 



|fc, r; k', r>) = e in |fc, r; k', r 1 ) « / 12 |fc, r; fc', r') (9) 

where /12 is the short range correlation factor, the so-called Jastrow factor 01 ■ F° r simplicity, we consider /12 = 
f(r 12), r\2 = f\ — r*2, and /(r) = 1 — (a + /?r) e~ ir where a, (3 and 7 are parameters. The important effect of the short 
range correlations is the replacement, in the expression for the ground-state energy, of the interaction matrix element 

(k,r;k',r'\Vi2\k, r;k',r') by (fc, r; fc', r'\V\2 +t\ + t2\k, r;k\r') — (fc, r; fc', r'\U + t2\k,r;k',r'), where t% is the kinetic 
energy operator of particle i. As argued by Moszkowski and Bethe |12|. it is expected that the true ground-state 
wave function of the nucleus containing correlations coincide with the independent particle, or Hartree-Fock wave 
function, for interparticle distances r > rh ea i , where r^eai ~ 1 fm is the so-called "healing distance" . This behavior 
is a consequence of the constraints imposed by the Pauli Principle. A natural consequence of having the correlations 
introduced by an unitary operator is the occurrence of a normalization constraint on /(r), 

(f(r) - 1) d 3 r = 0. (10) 



The correlated ground state energy becomes 

S = ^ ' k 2 dk [|fc|sin X (fc) + Mcos X (fc)] + ^P-pl + ^-p% 
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(27T) 



k 2 dk k' 2 dk' I [|fc|sinx(fc) + 2 M cos x(fc)l I(fc, fc') + \k\ sm X {k') J(fc,fc')} 



1 f kf 



{2-kY 



k dk k! dk' 



Mk,k') +cos X (fc)cosx(fc') J2 Bi(k,k') + sin X {k) sin x(fc') ^2 ^(M') 



i—a.uj %—g.uj 



(11) 



where Ai, Bi, Ci , / and J are exchange integrals. In the above equation for the energy density, the first 
term results from the kinetic contribution, the second and third terms come respectively from the a and lu di- 
rect contributions from the potential energy with correlations, the fourth from the exchange correlation contri- 
bution from the kinetic energy, and the last one from the a + lo exchange contributions from the potential en- 
ergy with correlations. The direct correlation contribution is zero due to 11U|) . The angular integrals are given 
by Ai(k,k') = B^k.k') = 2tt gf/in f£ dcos9 Fi(k,k' ', cos d), C l {k 1 k l ) = 2n g 2 /4n ^ cos0 dcosf? F^k, fc',cos0), 
I(k, k') — 2ir J? dcos8 C\{k, fc', cos#), and J(k, k') = 2tt J? cos 9 dcosO C\(k, fc', cos#), where 

Fi(k,k')= [ [f(r)V T (r)f(r)] e 1 ^'^ dr and C x {%, fc') = f(f 2 (r) - 1) e ^-*'K df (12) 



The baryon density and the scalar density are 



2 v f k '' - 2 v k 3 f 2v f k f - - 

PB = Wfl* dk = ^ 2 ~' Ps ^J^rJo cosx{k)dk - (13) 
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TABLE I: Parameters and ground state properties of nuclear matter at saturation density are given. These results were obtained 
with fixed: M = 939 MeV m CT = 550 MeV, m u = 783 MeV at k F0 = 1.3 far 1 with binding energy E B = e/p - M = -15.75 
MeV. We have used a density dependent parameter (HF+corr) 7 = 600 + 400&F /&fo MeV for the correlation. 
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K 


M*/M 


T/p B - M 


V d / P B 


Ve/ps 


T c /ps 








(MeV) 


(MeV) 


(MeV) 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


Hartree 


il. 079 


13.806 






540 


0.540 


8.11 


-23.86 






HF 


10.432 


12.223 






585 


0.515 


5.87 


-37.45 


15.83 




HF+corr 


4.4559 


2.6098 13.855 


-2252.448 


1000 


429. 


0.625 


15.95 


-73.12 


20.46 


19.96 



The couplings g ai g^, the meson masses, m a and m w and also three more parameters from the short range correlation 
function, a, (3 and 7 have to be fixed. The couplings are chosen so as to satisfy the ground state properties of the 
nuclear matter. We choose m a — 550 MeV and take m w = 783 MeV. The normalization condition 11U|) determines 
(3. We fix a by minimizing the energy. For the parameter 7 we consider a function which increases linearly with the 
Fermi momentum, of the form 7 = a\ + 02 /cf/^fo where a\ and a-z are free parameters. This is consistent with the 
idea that the healing distance decreases as kp increases. 
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Fig. 1 The correlation function f(r) 
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Fig. 2 The variational angle x(^) 



The correlation function f(r) is plotted in figure 1 as a function of the relative distance. The inclusion of correlations 
introduces an extra node in the ground-state wave-function contrary to what generally happens in non-relativistic 
calculations with a hard core. In this case the wave function has a wound. The quantities cos(%(/c)/2) and sin(x(fc)/2) 
are plotted in figure 2. They show how the interaction and correlations make the wave-function deviate from the free 
wave- function, represented by a dotted line. The correlated angle x(^) lies between the Hartrce-Fock(HF) and the 
free wavefunction angles. As a consequence, we will see that the correlated effective mass will not decrease so fast 
with density as the HF effective mass. 

In table [fl we have tabulated the parameters used in our calculation together with the compressibility K, the 
relative effective mass M* /M, the kinetic energy T/pb — M, the direct and exchange parts of the potential energy 
iy<i/ Pb and V e /pB respectively) with correlation and the correlation contribution to the kinetic energy T c /ps, all 
calculated at the saturation point. Notice that a HF calculation produces an EOS which is stiffer than the one 
obtained at the Hartree level. However, the inclusion of correlations gives a larger effective mass than both Hartree 
and HF calculations and a softer EOS. In fact, the contribution of direct and exchange correlation terms are of the 
same order of magnitude of the other terms in the energy per particle. Hence, they cannot be disregarded. 

We have computed the binding energies as function of the density for the Hartree, HF and HF+Corr and compared 
with the quark-meson-coupling model (QMC) [l3j and a non- linear Walecka model NL3 [1J|, as can be seen from fig. 
3. The inclusion of correlations make the equation of state (EOS) softer than Hartree or HF calculations. NL3 and 
QMC also provide softer EOS around nuclear matter saturation density but around two times saturation density, the 
EOS with correlations is softer than NL3. In figure 4 we plot the effective mass as a function of density. If correlations 
are included the effective mass does not decrease so fast with the increase of density as in a Hartree or, even worse, 
HF calculation. This explains the softer behavior of the EOS with correlations. 

We conclude referring that, although correlation effects in the Hartree and HF calculations are may be taken 
partially into account by a correct choice of the coupling constants, the explicit introduction of correlations has other 
effects such as softening the EOS. 
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Fig. 3 EoS for different parametrizations Fig. 4 Effective mass as a function of density 
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